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Photoproduction in Ultra-Peripheral Heavy Ion Collisions at STAR 
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We present recent STAR results on photoproduction in ultra-peripheral relativistic heavy ion collisions. In these 
collisions the impact parameter of the beam particles is larger than the sum of their nuclear radii, so that they 
interact via their long-range Coulomb fields. STAR has measured the production of p°(770) mesons in exclusive 
reactions as well as in processes with mutual nuclear excitation of the beam particles. We present results for the 
p° production cross section in Au-Au collisions at -^/s^ = 200 GeV for coherent as well as incoherent coupling. 
The dependence of the cross section on the p*^ rapidity is compared to theoretical models. We also studied the ratio 
of coherent p° to direct tt^tt" production as well as the p° helicity matrix elements and we observe interference 
effects in the p° production. In addition STAR has measured the production of p° mesons in d-Au collisions at 
y/Sj^j^ = 200 GeV and that of e^e~ -pairs in Au-Au at y^s^ 
1510 MeV/c^ in 7r^7r~7r^7r~ final states in Au-Au collisions at 



200 GeV. We also see an enhancement around 
= 200 GeV. 



On 
cn 

00 

o 

00 

o 



1. Introduction 

Heavy ions are surrounded by a cloud of vir- 
tual photons. The electromagnetic field of a rela- 
tivistic nucleus can be approximated by a flux of 
quasi-real virtual photons using the Weizsacker- 
Williams approach [1]. Since the number of pho- 
tons grows with the square of the nuclear charge, 
fast moving heavy ions can generate large photon 
fluxes. Beams of heavy ions can thus be used as 
photon sources or targets. Due to the long range 
of the electromagnetic interactions, they can be 
separated from the hadronic interactions by re- 
quiring impact parameters b larger than the sum 
of the nuclear radii Ra of the beam particles. 
These so-called Ultra-Peripheral heavy ion Colli- 
sions (UPCs) allow to study photonuclear effects 
as well as photon-photon interactions [2j . 

A typical photonuclear reaction is the produc- 
tion of vector mesons. In this process the vir- 
tual photon, that is emitted by the "spectator" 
nucleus, fluctuates into a virtual (7^-pair, which 
scatters elastically off the "target" nucleus, thus 
producing a real vector meson. The scattering 
can be described in terms of soft Pomeron ex- 
change. Figure [1^ illustrates this reaction taking 



as an example the most abundant process, the 
photonuclear production of p^(770) mesons. 

The cross section for vector meson production 
depends on the way the virtual qq-pdiii couples to 
the target nucleus. This is mainly determined by 
the transverse momentum pt of the produced me- 
son. For small transverse momenta of the order of 
Pt ^ 2h/RA the gg-pair couples coherently to the 
entire nucleus. This leads to large cross sections 
which depend on the nuclear form factor F{t) and 
which scale roughly with the atomic number A 
squared. For larger transverse momenta the qq- 
pair couples to the individual nucleons in the tar- 
get nucleus. This "incoherent" scattering has a 
smaller cross section that scales approximately 
with A modulo corrections for nuclear absorption 
of the meson. 

Due to the intense photon flux, it is possible 
for the vector meson production to be accom- 
panied by Coulomb excitation of the beam par- 
ticles. The excited ions mostly decay via the 
emission of neutrons which can be measured by 
the Zero Degree Calorimeters (ZDCs) [3 of the 
STAR setup. To lowest order events with mutual 
nuclear breakup are described by three-photon 
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Figure 1. (a) Exclusive photonuclear pro- 
duction in Au-Au ultra-peripheral collisions, 
(b) Photonuclear production accompanied by 
mutual Coulomb excitation of the ions. The 
p^ creation in 7*P fusion is adumbrated in form 
of a hatched circle. 



exchange (see Fig. [Ud): One photon to produce 
the vector meson and two photons to excite the 
two nuclei [4]. All three photon exchanges are 
in good approximation independent, so that the 
cross section for the production of a vector me- 
son V can be factorized: 

cr(AuAu ^ Au*Au* V) = 

J [1 - Phad(^)] Pv{b) PxnAh) PxnAb) (1) 

where Phad is the probability for hadronic in- 
teraction, Py the probability to produce a vec- 
tor meson V, and Pxn,i the probability that nu- 
cleus i emits X neutrons. Compared to exclusive 
vector meson production, reactions with mutual 
Coulomb excitation have smaller median impact 
parameters. 

In this paper we present results from the STAR 
experiment at the Relativistic Heavy Ion Col- 
lider (RHIC). Charged tracks were reconstructed 
using a large cylindrical Time Projection Cham- 
ber (TPC) [5J with 2 m radius and a length of 
4.2 m, that was operated in a 0.5 T solenoidal 
magnetic field. For tracks with pseudorapid- 
ity 1 7^ I < 1.2 and transverse momentum pt > 
100 MeV/c the tracking efficiency is better than 
85 %. For triggering two detector systems were 
used: the Central Trigger Barrel (CTB) [6j, which 
is an array of 240 scintillator slats surrounding 
the TPC that allows to trigger on charged par- 



ticle multiplicities in the event, and the two zero 
degree calorimeters, which are located at ±18 m 
from the interaction point and measure neutrons 
originating from nuclear breakup. 

2. Photonuclear Production 

STAR uses two kinds of triggers to select UPC 
events. The CTB-based "topology" trigger re- 
quires a low total charged particle multiplicity 
and divides the CTB detector into four azimuthal 
quadrants. Events are recorded when the left 
and right quadrants have coincident hits, select- 
ing roughly back- to-back pion pairs. The top and 
bottom quadrants serve as a cosmic ray veto. The 
"minimum bias" trigger selects events where both 
nuclei dissociated by requiring coincident energy 
deposits in both ZDCs. 

In the offline analysis the events were required 
to have two tracks of opposite charge that form 
a common (primary) vertex. In addition we de- 
manded a low overall track multiplicity in order 
to suppress backgrounds from beam-gas interac- 
tions, peripheral hadronic interactions, and pile- 
up events. The primary vertex was confined to a 
region close to the interaction diamond which re- 
duces contaminations from pile-up events, beam- 
gas interactions, and cosmic rays. The latter 
background component was decreased further by 
either excluding events around = in the 
topology data sample or by requiring neutron sig- 
nals in the ZDCs as in the minimum bias trig- 
ger. The contamination from peripheral hadronic 
interactions was in addition suppressed by de- 
manding the track pair to have low net transverse 
momentum. The backgrounds from two-photon 
e+e~ (cf. section [3j) and photonuclear oo produc- 
tion are negligible. 

The 7r+7r~ invariant mass distributions are fit 
with the function 



da 
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(2) 



Equation ([2|) consists of a relativistic Breit- 
Wigner with amplitude A, that describes the res- 
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onant tt+tt" production, and a mass-independent 
amplitude B for the non-resonant tt+tt" produc- 
tion, as well as a Soding interference term [7] 
of the two amplitudes. The second order poly- 
nomial, /bg, describes the combinatorial back- 
ground which was estimated from like sign pion 
pairs. 

2.1. Coherent Production in Au-Au 
Collisions 

Coherent photonuclear production is 

strongly enhanced at low p^. It can be separated 
from the incoherent production by requiring 
transverse momenta smaller than 150 MeV/c. 
Figure [2] shows the tt+tt" invariant mass dis- 
tribution of the minimum bias data set taken 
during the year 2002 run in Au-Au collisions 
at ^/s^ = 200 GeV [8J. The peak contains 
3 075 ± 128stat. candidates. 

Since in Eq. ([2]) the amplitudes A and B are 
free fit parameters, it is possible to extract the 
ratio of non-resonant to resonant tt+tt" 

production. In Au-Au collisions at ^/s^ = 
200 GeV the ratio was measured using the min- 
imum bias data to \B/A\ = 0.89 ± 0.08stat. ± 
0.09syst. GeV~^/^ [8]. As expected the ratio shows 
no dependence on the p^ rapidity nor on the 
p^ decay angles. The value is compatible with 
an earlier result of \B/A\ = 0.81 ± 0.08stat. ± 
0.20syst. GeV-^/2 from the year 2000 run at 
^ys^ = 130 GeV [9j. Both values are also in 
agreement with the ZEUS measurements [10 . 

The luminosity of the minimum bias sample 
was measured based on the known total hadronic 
production cross section of 7.2 barn. With the 
value of Lminbias = 461 mb~^ the cross sec- 
tion for coherent p^ production accompanied by 
mutual nuclear excitation was estimated to be 
14.5 ± 0.7stat. i l-9syst. mb~^ within the exper- 
imental acceptance of |?/p| < 1 

There are at least three models for p^ produc- 
tion in ultra-peripheral collisions: The model of 
Klein and Nystrand (KN) pjj uses the Vector 
Dominance Model (VDM) to describe the virtual 
photon and a classical mechanical approach for 
the scattering on the target nucleus, based on re- 
sults from jp p^p experiments. The Frank- 
furt, Strikman, and Zhalov (FSZ) model [12j em- 




Inv. Mass, GeV/c 

Figure 2. tt+tt" invariant mass distribution for 
coherently produced p^ from the minimum bias 
data set [8]. The solid line is the fit result using 
Eq. ([2]), which has four components: The rela- 
tivistic Breit-Wigner of the resonant tt+tt" pro- 
duction (dashed line), the mass-independent term 
for the non-resonant tt+tt" production (dash- 
dotted line), the interference of the two (dotted 
line), and the combinatorial background (hatched 



ploys a generalized VDM for the virtual pho- 
ton and a QCD Gribov-Glauber approach for 
the scattering. The model of Goncalves and 
Machado (GM) ^ is based on a QCD color 
dipole approach that takes into account nuclear 
effects and parton saturation phenomena. Only 
the KN model gives predictions for the rapidity 
dependence of the coherent p^ production cross 
section accompanied by mutual Coulomb excita- 
tion, which are compared to the data in Fig. [3l 

The total coherent p^ production cross section 
was estimated using the 13 054 ± 124stat. can- 
didates in the data set taken with the topology 
trigger [8j, which does not apply any restrictions 
on the ZDC amplitudes. Although this data sam- 
ple could not be used for an absolute cross sec- 
tion measurement, because the trigger efficiency 
is only poorly known, it was possible to extract 
coherent cross section ratios for the different nu- 
clear excitation states. These ratios were used to 
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Figure 3. Differential cross section da/dyp for 
coherent production accompanied by mutual 
Coulomb excitation as a function of rapidity 
extracted from the minimum bias data The 
solid line shows the KN model calculation [11]. 



scale the coherent mutual-excitation cross section 
values from the minimum bias sample. Figure [H 
compares the obtained ^p-dependence of the co- 
herent p^ production cross section to the three 
models. Unfortunately the limited acceptance in 
yp does not allow to distinguish between different 
shapes of the distributions. 

STAR has also measured the spin structure of 
the coherent p^ production amplitudes in terms 
of spin density matrix elements [8j. The matrix 
elements were extracted from the minimum bias 
data by fitting the p^ decay angular distribution 
in the p^ rest frame to the theoretical prediction 
for unpolarized photoproduction [M] 



1 dV 



a dcos^ 



3 

47r 



1 



-(3^00 - 1) cos^ - V2D\t[r^^] sin 2(9 cose 



r?ii sin^6> cos 2^ 



(3) 



where is the polar angle between the beam di- 
rection and the direction of the 7r+ and (j) the 
azimuthal angle between the decay and the pro- 
duction plane. Since STAR does not measure 



Figure 4. Comparison of the measured cross sec- 
tion for coherent p^ production with theoretical 
predictions [8j. The vertical line at each point 
shows the statistical error. The shaded area dis- 
plays the sum of statistical and systematical er- 
rors. The dashed line represents the KN [11], the 
dash-dotted line the FSZ [12], and the dotted one 
the CM model [13]. 



the scattered beam particles, only an approxi- 
mate production plane with respect to the beam 
direction could be reconstructed. Tqq represents 
the probability that the p^ is produced with he- 
licity 0, ri^_i the interference between helicity 
non-flip and double-flip. Due to the ambiguity 
of photon source and target, 5^e[rJo], which is 
related to the interference between helicity non- 
flip and single-flip, could not be measured. Ta- 
ble [T] summarizes the results. The measured spin 
density elements are small, indicating 5-channel 
helicity conservation. They are also compatible 
with ZEUS measurements on p^ photoproduction 
from proton targets [10]. 

2.2. Interference Effects in Coherent 
pO Production in Au-Au Collisions 

In ultra-peripheral collisions the impact param- 
eter b is typically much larger than the sum of 
the nuclear radii of the beam particles. The 
scattering of the virtual qq-pdiii on the target 
nucleus, on the other hand, involves the short- 
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Table 1 

Spin density matrix elements for coherent production accompanied by mutual Coulomb excitation in 
Au-Au collisions at -^/s^ = 200 GeV compared to ZEUS results from production in collisions [lOj. 



Parameter 






STAR 










ZEUS 








-0.03 


± 


O.OSstat. 


± 


0.06syst. 


0.01 


± 


O.Olstat 


± 


0.02syst. 














0.01 


± 


O.Olstat 


± 


O.Olsyst. 




-0.01 


± 


O.OSstat. 


± 


O.OSsyst. 


-0.01 


± 


O.Olstat 


± 


O.Olsyst. 



ranged (0(1 fm)) strong interaction, so that the 
production of the p^ mesons occurs in or very near 
to the target nucleus. Furthermore photon source 
and target are indistinguishable, so that either 
nucleus 1 emits a photon which produces a p^ at 
nucleus 2, or vice versa. The system thus behaves 
like a two-source interferometer with separation 6, 
where the interference between the two possible 
amplitudes creates an entangled final state tt+tt" 
wave function. Since the p^ has negative parity, 
the two amplitudes have different sign, so that 
the cross section can be written as [15] 



a{b) = A{b,yp)-A{b,-yp) 



_ 2 



(4) 



where A{b^ yp) and A{b, —yp) are the p^ produc- 
tion amplitudes for the two photon directions. 
The observed pj. spectrum is obtained by inte- 
grating Eq. dl]) over 6, which is not measured in 
the experiment. For yp = the two amplitudes 
are equal so that 



cr{b) = cro{b) 1 — cos{p^ • b) 



(5) 



based on the KN model [llj. A is the overall 
normalization, B the slope parameter for coher- 
ent p^ production, and c measures the strength of 
the interference: c = corresponds to no interfer- 
ence, whereas c = 1 is the expected interference 
in the KN model. 
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Figure 5. Efficiency-corrected t-spectra for (a) 
the minimum bias data set with \yp\ < 0.5 and 
for (b) the topology data set with 0.05 < \yp\ < 
0.5. The minimum bias data exhibit a wider dip 
at low because the median impact parameter is 
smaller. The solid line shows a fit to Eq. (|6]). 



with ao{b) being the cross section without inter- 
ference. From Eq. ([5]) it is clear that the cross 
section is suppressed for p^ < 2h/ (6), where (b) 
is the mean impact parameter. 

The t-spectra, with t ~ (p^) , are parametrized 

by 



:^=Ae-^'[l + c{R{t)-l)] 



where = [dN / dt]^^ / [dN / dtj^^^, is the 
ratio of the ^-distributions with and without inter- 
ference obtained from Monte-Carlo simulations 



(6) 



iMC 



Due to the neutron tagging of mutual Coulomb 
excitation in the minimum bias trigger, the me- 
dian impact parameter of b ~ 18 fm in this 
data set is smaller than the b ^ 46 fm in the 
topology data. Therefore the dip in the mini- 
mum bias t-spectra should extend to larger val- 
ues of t. Figure [5] shows the efficiency-corrected 
^-distributions for the run 2002 topology and min- 
imum bias data sets from Au-Au collisions at 
y/s^ = 200 GeV in the central rapidity region 
\yp\ < 0.5 together with a fit to Eq. (|6|). To avoid 
contamination of the topology sample by cosmic 
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rays, in addition the region \yp\ < 0.05 was ex- 
cluded. 

The minimum bias data exhibit an interference 
strength of c = 0.92 ± O.OTgtat. close to the ex- 
pected value. The value from the topology data 
set is lower with c = 0.73 ± 0.1 Ogtato but could be 
affected by the imperfect topology trigger simu- 
lation. 

2.3. Incoherent Production in Au-Au 
Collisions 

The incoherent production, where the 
couples to the individual nucleons, can be mea- 
sured by extending the analyzed range of the 
p^ transverse momentum to 550 MeV/c. If the re- 
gion of t < 0.002 (GeV/c)^, where interference ef- 
fects reduce the cross section (see subsection 12. 2p . 
is excluded, the ^-distribution can be described by 
a double exponential function 

^ = A,ohe-^-'>* + A„ce-«-* (7) 

Bcoh is the slope parameter of the exponen- 
tial that dominates the \ow-p^ region and corre- 
sponds to the coherent production. Accordingly 
Bine is the slope parameter of the incoherent pro- 
duction that extends to larger p^. 

Figure [6] shows a fit of Eq. ([7j) to the year 
2002 minimum bias data from Au-Au at ^/s^ = 
200 GeV that gives ^coh = 388±24stat. (GeV/c)-^ 
and Bine = 8.8 ± l.Ostat. (GeV/c)-^ '[8|. By in- 
tegrating the exponentials the cross section ratio 
of incoherent to coherent p^ production accompa- 
nied by mutual Coulomb excitation was measured 

to 29 ± 3stat. i 8syst. %• 

2.4. Production in d-Au Collisions 

In d-Au ultra-peripheral collisions the virtual 
photon is predominantly emitted by the gold nu- 
cleus, due to its larger charge. There are two 
possible p^ production processes: coherent scat- 
tering on the entire deuteron, j^d ^ p^d^ and 
"incoherent" scattering on the individual nucle- 
ons of the deuteron in the course of which the 
deuteron dissociates, 7*6/ p^ pn. 

During the year 2003 run STAR took data of 
d-An collisions at ^/s^ = 200 GeV using two 
UPC triggers. The first trigger was a topology 




0.05 0.1 0.15 0.2 0.25 0.3 

t, (GeV/c)2 



Figure 6. Dependence of the p^ production cross 
section on the squared momentum transfer t for 
the minimum bias data [8j. The solid line shows 
the result of a fit to the double exponential Eq. ([7j) 
with the steep line at low t corresponding to co- 
herent production. The line at larger t describes 
the incoherent p^ production. 



trigger which selected mainly coherent interac- 
tions. The second trigger required, in addition 
to the topology condition, a neutron signal from 
the deuteron breakup in the corresponding ZDC. 
This data sample contains incoherent scattering 
events. 

Figure [71 shows the t-spectrum, with t ^ 
{p^) , for incoherent scattering accompanied by 
deuteron dissociation. A fit with an exponen- 
tial gives a slope parameter of B = 9.1 ± 
0-9stat. (GeV/c)~^ which is related to the nucleon 
form factor. This value is compatible with the 
slope parameter of Bine = 8. 8± l.Ostat. (GeV/c)~^ 
measured in incoherent p^ production accompa- 
nied by mutual Coulomb excitation in Au-Au 
collisions at ^/s^ = 200 GeV (cf. subsec- 
tion \2.3\ [8J) and also agrees with the value of 
B = 10.9±0.3stat.^o:5syst. (GeV/c)-^ measured 
at ZEUS [lOj. The downturn of the distribution 
at low t, which is also seen in low energy jd ex- 
periments [1^ , is due to the fact that there is not 
enough energy to dissociate the deuteron. 
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Figure 7. t-distribution for incoherent produc- 
tion accompanied by deuteron dissociation. The 
sohd hue shows an exponential fit. 



3. e^e -Pair Production in Photon- 
Photon Interactions 

During the year 2002 run STAR also measured 
the production of e+e~ -pairs in photon-photon 
interactions in ultra-peripheral Au-Au collisions 

V^NN — 200 GeV [17j. Since the electromag- 
netic field is strong with Za ~ 0.6, one expects to 
see higher-order QED effects. In order to enrich 
collisions at small impact parameters, where the 
fields are strongest and thus higher-order correc- 
tions should be most pronounced, a ZDC trigger 
was employed, requiring mutual nuclear dissoci- 
ation and thus limiting the impact parameter to 
2i^A < ^ < 30 fm. 

In the offline analysis events were selected that 
had two reconstructed tracks of opposite charge 
emerging from a common (primary) vertex along 
the beamline. Up to two additional tracks, not 
associated to the primary vertex, were allowed in 
order to account for random backgrounds. Most 
of the are produced in forward direction with 
low transverse momentum, so that only a small 
fraction of the pairs reached the central detectors, 
even though the solenoid field was reduced to 
0.25 T, half of its nominal value. The tracks were 
required to have pseudorapidities of |7^| < 1.15 
and transverse momenta of 65 < pt < 130 MeV/c. 
The latter cut ensured that the identification 
via the specific energy loss in the TPC gas had 



only minimal contaminations from other parti- 
cle species. Additional cuts on the pair trans- 
verse momentum, ^ < 100 MeV/c, and the 
pair rapidity, |^e+e-| < 1.15, gave a final sample 
of 52 events in the pair invariant mass region of 
140 < Me+e- < 265 MeV/c^ 

The measured cross section in the selected kine- 
matic range is 1.6±0.2stat.i0.3syst. mb. Figure [8^ 
shows the cross section as a function of the trans- 
verse momentum of the e+e~-pair. Fig. [SJd the 
pair invariant mass distribution. The data are 
compared to a a recent QED calculation of Baltz 
that includes a realistic phenomenological treat- 
ment of the nuclear Coulomb excitation [18j. The 
lowest order QED (dashed line in Fig. [8]) over- 
predicts the data. The calculated cross section 
of 2.34 mb is nearly 2 standard deviations above 
the data. The corresponding higher-order QED 
calculation gives a cross cestion value of 1.67 mb, 
compatible with the measurement, and also de- 
scribes the kinematic distributions well. 




Pair Transverse Momentum (MeV/c) Pair Invariant Mass M,, (MeV) 



a) b) 

Figure 8. (a) Pair transverse momentum spec- 
trum, (b) e+e~ invariant mass distribution. The 
dashed line shows the lowest order, the solid 
line the higher-order QED calculation according 
to [18J. 



4. Coherent Photonuclear tt+tt tt+tt 
Production 

During the year 2004 run STAR took data in 
Au-Au collisions at y/s^^^ = 200 GeV using a 
"multi-prong" trigger which combined the UPC 
neutron tag in form of coincident ZDC signals 
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with a low hit multiphcity in the CTB. In addi- 
tion a veto from the large-tile Beam-Beam Coun- 
ters (BBCs) [19], which cover 2.1 < |7^| < 3.6, was 
applied. 

Events were selected that have four charged 
tracks coming from one common (primary) ver- 
tex. These four-prongs were required to be neu- 
tral and to have a net transverse momentum be- 
low 150 MeV/c, thus selecting coherent produc- 
tion. A low number of additional tracks, unasso- 
ciated to the primary vertex, were allowed in the 
events to account for backgrounds. 




Figure 9. (a) Transverse momentum spectrum 
of coherently produced neutral four-prongs. The 
points are the data, the shaded histogram depicts 
the expected distribution from simulation. The 
hatched histogram shows the background esti- 
mated from charged four-prongs, (b) 7r+7r~7r+7r~ 
invariant mass distribution of coherently pro- 
duced four-prongs. A broad peak indicates res- 
onant 7r+7r~7r+7r~ production. 



The transverse momentum distribution of the 
neutral four-prongs in Fig. [9^ shows the expected 
peak at low pt, characteristic for coherent pro- 
duction. The 7r+7r~7r+7r~ invariant mass spec- 
trum in Fig. [Hb shows a broad peak of approx- 
imately 100 events around 1510 MeV/c^ which 
could be caused by the production of excited 
states like the p°(1450) and/or the p°(1700). 
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